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O 

f"^ ' The roughness of vapor-deposited thin films can display a nonmonotonic dependence on film thick- 

ness, if the smoothening of the small-scale features of the substrate dominates over growth-induced 
roughening in the early stage of evolution. We present a detailed analysis of this phenomenon in the 
» ! ' framework of the continuum theory of unstable homoepitaxy. Using the spherical approximation of 

f> |, phase ordering kinetics, the effect of nonlinearities and noise can be treated explicitly. The substrate 

■ roughness is characterized by the dimensionless parameter Q — Wo / (koa 2 ) , where Wo denotes the 

roughness amplitude, fco is the small scale cutoff wavenumber of the roughness spectrum, and a is 
the lattice constant. Depending on Q, the diffusion length In and the Ehrlich-Schwoebel length 
Ies, five regimes are identified in which the position of the roughness minimum is determined by 
different physical mechanisms. The analytic estimates are compared by numerical simulations of 
i <~| the full nonlinear evolution equation. 
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I. INTRODUCTION 



The morphology of thin film surfaces has a decisive influence on many film properties. The control of growth- 
induced surface roughness is therefore a central concern in thin film science and technology. Two types of roughening 
mechanisms have been extensively studied in recent years . The term kinetic roughening is commonly used to refer 
to a stochastic mechanism, in which fluctuations in the deposition flux interact with thermal smoothening to generate a 
scale-invariant, rough morphology. This theoretically appealing but empirically rather elusive phenomenon |3j is often 
superseded by a second, deterministic mechanism, a growth instability associated with reduced interlayer transport 
and slope-dependent mass currents along the surface y||J| ■ The hallmark of unstable growth is a morphology of more 
or less regular mounds with a clearly developed characteristic length scale. While in practice the distinction between 
the two types of roughening mechanisms may not always be so clear-cut || , they are very different conceptually. 

In addition to the growth-induced roughness, clearly also the roughness of the substrate affects the film morphology. 
Since the growing film covers up the small-scale details of the substrate modulations, the substrate contribution to 
;> the roughness is expected to decrease with increasing film thickness, while the growth-induced roughness component 
increases. Under suitable conditions this leads to the somewhat counterintuitive possibility of a minimum, of the total 
surface roughness at a nonzero film thickness. This phenomenon has been observed in several growth experiments 
[ , and a theoretical description has been worked out on the level of linear continuum theories of kinetic roughening 
[ HJ and unstable growth . 

In Ref. ]TT[ | a quantitative comparison with the experiments of Gyure et al. |)| was attempted, which indicated an 
important influence of nonlinearities. This motivated the present study, in which the nonlinear term in the growth 
equation is treated explicitly using the spherical approximation of phase ordering kinetics p^ , ^3| . We find that the 
interplay of instability, nonlinearity and noise gives rise to a rather complex behavior, in which the position of the 
roughness minimum can be determined by several distinct physical mechanisms. For a quick overview of the different 
I 1 regimes we refer the curious reader to Table I. 

The paper is organized as follows. In the next section we introduce the standard continuum equation for unstable 
homoepitaxial growth Jl3|,[l4j and describe the strategy for its analytical solution. Section |n|is devoted to the 
roughness evolution in the absence of noise. We first recapitulate the linear analysis of Ref. flll|, then provide a 
detailed analysis of the relevance of the nonlinearity and the nonlinear behavior, and finally discuss the influence of 
correlated initial roughness. The effects of noise are analyzed in Section |v|. In Section [v| we compare the analytic 
estimates to a numerical evaluation of the spherical approximation, as well as to numerical simulations of the full 
nonlinear growth equation, finding good agreement in all cases. Finally, some conclusions are formulated in Section 
" |VT. 



II. THE CONTINUUM EQUATION 

The evolution of a surface growing under typical Molecular Beam Epitaxy (MBE) conditions is described by an 
equation of the form [Q 

d t H + V • J = a ± F + n, (1) 
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where H is the height, J is the surface current, a± is the monolayer thickness, F is the average value of the deposition 
flux and r\ is a noise term describing fluctuations in the flux (shot noise) and in the diffusion of adatoms on the surface. 
The constant term a±F can be eliminated passing to the frame of reference H = a±Ft + h moving with the average 
height. The noise has zero average and correlations 

Mx, tMx',0) = (2n)- 2 6(t - t')(R s - i? D V 2 )<5(x - x'), (2) 

with the amplitudes R$ and Rd representing the effect of shot and diffusion noise, respectively. 
The surface current is the sum of two contributions 



J = kV(V 2 /i) + f[(Vh) 2 ]Vh. (3) 

The first tends to smoothen the surface and has the form of a capillarity term jl5| ] , even though it may be dominated 
by nonequilibrium effects such as nucleation [ h4 | |lq ]. The second term models a growth- induced surface current, whose 
existence is often (but not necessarily jlTHl 9ff) caused by the presence of an Ehrlich-Schwoebel barrier for interlayer 
diffusion. We assume in-plane isotropy of the current (for a discussion of origins and consequences of anisotropy see 
IP^POI ). For the function /(a) we use the form valid in the weak barrier limit || 

= (l + Z DV ^/a±)(l + W^/ol)' (4) 

where lp is the diffusion length and Ies the Ehrlich-Schwoebel length. These length scales are related to the in-layer 
hopping rate D, the inter-layer hopping rate D' and the deposition flux F through 

l D » (D/Fya {] , l ES = (D/D 1 - ljoy . (5) 

Here a\\ denotes the in-layer lattice constant and the exponent 7 depends on the size of the critical cluster for two- 
dimensional nucleation plj . Conditions of weak and strong step edge barriers can be distinguished according to 
whether Ies 3> Id (strong barriers) or Ies <C Id (weak barriers) P,|^,[l6f . Here we focus on the latter case, in which 
a continuum description is most likely to be valid. 

The coefficients a and k in Eqs.(|[Q) are related to microscopic parameters by a w FIesId/2, « ~ Fl D P,|5|, [r^ , ^6| . 
We will for simplicity assume that the equality sign holds in these formulas; however numerical factors are not precisely 
known and the equalities should be intended instead only in order of magnitude. Also the amplitude of the noise 
terms is connected to microscopic parameters through Re — ^-Rg = if^Fa^a 2 . 

By inserting the expression ([}]) of the current in the equation for the height and neglecting noise one obtains 

«9 t /i = - K (V 2 ) 2 /i-V[/((V/i) 2 )V/i]. (6) 

This strongly nonlinear equation is reminiscent of the Cahn-Hilliard equation for phase-ordering in systems with 
conserved order parameter fl^| . A widely used method for investigating this kind of nonlinear evolution is the large- 
limit, or spherical approximation. In the present context it consists in replacing the argument of the nonlinear 
current / with its average value a(t) = ((Wh) 2 ). In this way, Eq. (||) is effectively linearized. It is then possible to 
write down a closed form linear equation for the structure factor S(k,t) = (h(k. ) t)h(—'k,t)) 

d t S(k, t) = -2 [ K k 4 - f[a{t)]k 2 ] S(k, t), (7) 

where a(t) must be determined self-consistently 

a{t) = 2ir dkk 3 S(k,t). (8) 
Jo 

The solution of this pair of coupled equations has already been derived for long times by Rost and Krug [ p^[ . 
Here we concentrate on the short time behavior, i.e. all what happens before the instability sets in. In this time 
range we expect the large- A approximation to give a fairly accurate description of the nonlinear behavior ]22| , since 
correlations are still small in range and amplitude. In particular we will be interested in the time evolution of the 
surface roughness 

W 2 (t) = 2ir dkkS(k,t). (9) 
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We usually assume as initial condition So{k) = S(k = 0, t) a white spectrum with an upper cutoff fco = tt/Iq 

which implies W 2 (0) = W§ and a(0) = WqIiq/2. The dimensionless number 

Q = Wo/(k a x a l{ ) (11) 
will turn out to provide a useful measure for the strength of the initial roughness; note that it involves b oth th e 



amplitude (Wq) and the small scale cutoff. Other types of initial roughness spectra will be treated in Section [II C 



III. SOLUTION IN THE DETERMINISTIC CASE 



Equation ([?]) can be formally integrated 

S* det (fc, t) = 2 f 2^ exp [~2nk 4 t + 2k 2 b(t)] 



where 



b(t)= / dsf[a(s)}. 
Jo 



By defining k^(t) = b(t)/(2nt) one can rewrite 



W 2 {t) 



Wo 



dk k exp < 2tnk 



In lA ^ [Erfjfcgygjg + Erf((fc 2 - kpy/M)] 

W x j - eM2 K k m t) 



(12) 
(13) 

(14) 
(15) 



with Erf(s) = (2/y/n) J* exp(— t 2 )dt. The wavenumber k m is the position of the structure function peak (when it is 

a real number, otherwise the peak is for k = 0). Within linear theory its value is constant, k m — ki = y // a/2K, while 
in general k rn is a function of time. 

The derivative of W 2 with respect to t at t = is 



dW 2 {t) 



At 



'-W 2 K k 2 (3~k 2 -k 2 ), 



(16) 



where k = fc ro (0) = v //[a(0)]/2« = y/f(W§ k 2 /2)/2n. Hence if k/k > l/%/3 the roughness grows from the beginning 
and there is no minimum in the behavior of W 2 (t): The instability is immediately at work. Notice however that the 
condition for the existence of the minimum involves k and not fc;. When the initial roughness is large k <C fc; and it 
may occur that k <C fco "C kf. In such a case a minimum occurs even if the linear theory does not predict it. This fact, 
together with the observation that h/ko does not depend on Wo, implies that if fco is sufficiently large a minimum 
exists even for very small amplitude of the initial fluctuations. If fco is not large, only a strong initial roughness can 
originate a non monotonic behavior of the width. 

We will assume in the following that a minimum exists. To study the detailed behavior of the system one should 
in principle consider Eq. (|l^) and (|l3|) simultaneously. However, expanding the expression of W 2 (t) for small k m /ko 



w< it) - ( £ 



TrErf^V^) 



2Jnt 



W 2 [l-exp(-2«fc^)] ( J?L 



fco 



(17) 



one finds that after a transient time to — l/(2KfcQ) the width starts decreasing as l/v^i and it does so until 1/vKt — 
fc^. During this time interval the width decreases in time and depends only on k and fco, not on k m (hence not on the 
form of the current, not even its linear expansion). The system is effectively described by an evolution equation (0) 
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where only the relaxational term proportional to k matters |ll],^3j . This fact is crucial for all the following calculations: 
The structure factor is known 

S det {k,t) = So cxp(-2 K fc 4 i) (18) 

and can be used to compute a(t) and k m (t), which change nonlinearly in time, but do not affect significantly Sd e t(t)- 
This situation persists up to a time such that other terms in the expansion (|l^) become large. Since the other terms 
grow with t, it is around this time that W 2 (t) reaches a minimum. The only role played by the current / is to 
determine the time evolution of k m (t) and hence when the decay of the initial fluctuations ends: But / does not affect 
the way W 2 (t) decreases. 

The form rtlq) of the structure factor may be interpreted as that of a system which is coarsening with a typical 
correlation length growing as L(t) ~ (nt) 1 / 4 . The initial condition creates "domains" of size k^ 1 , much smaller than 
the length of the instability k^ . the system evolves by reducing the amplitude of fluctuations and increasing the 
correlation length. This explains why the initial decrease of W 2 is seen only for small k/ko- This coarsening process 
continues until 

L^k~\ (19) 

^From this time on the evolution proceeds by amplifying fluctuations of scale L and the instability sets in. Notice 
that while in the linear case k m = fc; is constant in time, when nonlinearities are taken into account k m grows in time 
but remains always smaller than ki, because /(a) < a for all a. Therefore the time where the minimum occurs in the 
linear theory is a lower bound for the same quantity in the nonlinear case. 

We now compute in detail how the position of the minimum depends on the amplitude of the initial fluctuations. 
For reference it is useful to summarize first the results of the linear theory |ll| . 

A. Linear theory 

The assumption of linearity for the current implies a(t) = 0. Hence f[a(t)] = f[0] = a, b(t) = at and k^ = kf = 
a /2k. Then the temporal evolution of the width is fully specified for all times by Eq. (|l5|). Letting k — > oo such a 
formula can be cast as 

W 2 (t) = ^^\(t/nl (20) 

where 

ti = An/a 2 (21) 

is the inverse amplification rate of the maximally unstable fluctuations and the scaling function is 

$(x) = e 2x ^j2-K/x[\ + Erf(V2aO]. (22) 

This formula is valid only for times greater than to — Xjlnk 4 ~ (fc;//co) 4 T; <C V- 
The width attains a minimum at a time 

4^ « 0.18 ti, (23) 

where it has been reduced by a factor 

W 2 (t l mm )/W 2 « 3.42 (h/ko) 2 . (24) 

This minimum marks the transition between the initial power-law decrease and the eventual exponential increase due 
to the linear instability. 
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B. Nonlinear theory 



The initial condition implies that a(0) = W^k^/I. Then, for very small t 

Ht) = f(W%kl/2)t = at. (25) 
By comparing a with a one recovers the condition (21) of Ref. [p"lf for the irrelevance of the nonlinearity 

W k l D 



< 1. 



(26) 



However, this condition turns out to be too restrictive. The reason is that it assesses the relevance of the nonlinearity 
from its importance in the expression of the unstable current at the initial time. But, as discussed above, / does not 
play any role in the initial evolution of the width. The relevance of the nonlinearity must instead be established from 
its influence on the position of the minimum of W 2 {t), i.e. for long times. By that time the initial fluctuations have 
already been reduced significantly. 

Expanding a(t) as a function of k m /ko 



a(t) 



wl_ 

Anklt 



[1 - cxp(-2Kfcot)] 



/2tt 



-4exp(-2Kfcpt)fco + \ — Erf(fc„V2Kt) 

" Kt 



km 

k 



()[ 'km 

k 



(27) 
(28) 



one can see that a(t) is constant only for times of the order of to = l/(2«fcg), indicating that Eq. ( |2q ) soon loses 
validity. 

For longer times (up to t ~ 1 / [Kk m (t) 4 )) , 



a(t) 



Wp 2 
Anklt 



Inserting Eq. (E9T) in the expression (O) we can compute the long time behavior of b(t) 



b(t) = 



ds 



(1 



lpW 
2a ± ko^/K,s 

a 



)(i 



IesWq ' 
2a±ko^ ks ' 



= ds ■=. . 

(l+v^)(l + vw) 



where we have introduced two timescales 



td 



IpWpk; 
2a±ko 



n = QV ( zr- 



a 

81d 



n 



(29) 

(30) 
(31) 

(32) 



and 



1~ES 



lEsWpkf 
2a±ko 



Tl = qV 



a ii 
8^ 



(33) 



In the right equalities we have introduced the quantity y = Ies/Id and Q is defined in (|l l|) . It is important to 
stress that since we are considering weak Ehrlich-Schwoebel barriers y <C 1 and hence tes "C td. Moreover, realistic 
values of the diffusion length are such that Id/o,\\ 3> 1. Depending on the initial roughness via the parameter Q, the 
timescales tes and td will be larger or smaller than the initial one, r;, giving rise to different scenarios. 



1. Irrelevant nonlinearity (intermediate initial roughness) 

Consider first the case of a fairly small initial roughness, i. e. Qya»/8lr) <C 1, so that tes *C td <C t\. Then for 
to -C t -c tes one can neglect the constant term in the denominator of Eq. (|3~i|) 
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as at 2 



For tes <t< instead 



y/ T D T ES 2y/T D T E S ^ ^ 



, s f* as 1 ' 2 2at 3 / 2 , s 

b(t) « / ds —= = —=, (35) 

while for td <C £ <C r; 

6(t) « at. (36) 

6(i) undergoes several changes during the time evolution passing through two intermediate behaviors. However it 
is easy to see that these variations in the form of b(t) are too short lived to affect the time evolution of W 2 (or a(t)). 
The minimum width is reached when the effect of the nonlinearity is already lost, and is well described by linear 
theory. The condition for the irrelevance of the nonlinearity is therefore that To <C n , that is 

^ (£)'«'■ 

Comparing with Eq. ( p6| ) it is clear that the relevance of the nonlinearity is strongly reduced when fc; <C kg. Notice 
moreover that the role of k is opposite compared to the condition ( p6] ) : An initial substrate rough down to very small 
length scales (large fco) makes the nonlinearity less relevant. 



2. Relevant nonlinearity (large initial roughness) 

Let us assume instead tes C^K td, i.e. 8lE>/(ya\\) <C Q <C 8lD/(y 2 a\\)- One has for tes <t <C id 

b(t) = -^t 3 / 2 « at. (38) 
3v r £> 



With this expression of b(t) the value of k m is much smaller than the linear value ki and the condition ( |19| ) for the 
minimum is attained on time scales larger than 77. To estimate t m i n more precisely the form of W 2 (t) can be written 
in the scaling form 

W 2 (t) = ^(^-)\(t K k m (t)% (39) 



4 \ k 

where $ is defined in (|2^). Now the minimum value of W 2 is not reached where $'(x) = 0, because k m depends on 
t. The condition for the minimum of W 2 is instead 

x$'(x)(2b't -b) + <f>(x)(b't - b) = (40) 

Using the expression ( |38| ) for b(t) one gets 

x®(x) 1 



<$>{x) 4 

whose solution is x mm = t° in /T(t^ m ) « 0.06 yielding 



C n «0.4(r^)^ = .4Q^|l) 



(41) 



n (42) 



Hence t® in is larger than 17 but smaller than td. The value of the width at the minimum is 

1/2 

W(t° in )/W « 1.3 (^J (43) 
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If instead Q 3> 8lD/(y 2 a\\), then 77 <C tes Ctd. By means of an analogous procedure, one finds that the minimum 
occurs for x m i n ~ 0.03 yielding 

« 0.5 (t,^) 1 / 8 = 0.5 Q^y 2 (44) 
The value of the width at the minimum is 



C. Correlated initial conditions 



The previous results can be easily extended to the case of a substrate with correlated roughness, i.e. with 

(A for k<k* 

S (k) = I A{k*/k) e for k* <k<k Q (46) 
[ for k > ko, 

with 6 > 0, and k /k* > 1. 

In this case the condition for an initial decrease of the roughness is, in the limit fc /fc* — ► oo, 

(*X<twh f ° r 6<A (47) 
\k ) \0 for 6»>4 k ' 

Hence for 6 > 4 the width can only increase monotonically from the beginning. 



1. Linear theory 



For < 2 one can safely take k* — > and find, for t 3> l/(2/cfco) and to second order in fc; 



I4^ 2 (i) =^fc* e (2Ki) e / 4 



(8Kt)V2 



2-6> 



fcfr [ i 



Estimating the position of the minimum from the time when the second term equals the first, one obtains 

-i 2 



t: 



l (0) 



1 



8nkf 



r(i-I) 



while the minimum width reached is 
W\t l min {6)) 



9/2 



rll -4 



1-61/2 



2-6 



(48) 



(49) 



(50) 



Both t l min (8) and W 2 (t l min (0)) are growing functions of 9: The minimum is delayed and made shallower by the 
presence of correlations in the roughness. This occurs because the roughness is concentrated on large length scales 
and the damping of small scale fluctuations provided by the rclaxational dynamics is less effective in the reduction of 
the surface width. Such effect is most evident when 6 approaches 2: t l min {6) diverges, while W 2 (t l min (0)) /Wq goes to 
1, since for 9 = 2 all the roughness is concentrated on the macroscopic length scale (A;*) -1 . 

For 2 < 9 < 4, the origin of the minimum is different. In this case one can take fco — > oo, and find for t <C 1/ (2«;fc* 4 ) 
and small fc; 



W 2 {t) = W$ + 2irAk* e/4 (2Kt) 



6/4 



1 



+ 



2k? 



(6 - 2)(2/si)V2 4-0' 



(51) 
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Differently from Eq. (fig), here t has, in the first contribution to the second term, a positive (and small) exponent but 
a negative prefactor. Therefore in this case the initial decrease of the roughness is much weaker and this is reflected 
by the minimum width, which is close to Wq- The time when the minimum is reached is 

and vanishes, as expected, in the limit 9 — > 4. 



2. Nonlinear theory 

With correlated initial conditions the evolution of the average square slope, considering S(k, t) = So(k) exp(— 2/-cfc 4 i) 
and k* -> 0, is, for t > 1/(2kA$) 



o(t) = 



irAk* e r{l - 9/4) 



(53) 



Hence a(t) is for all 9 < 4 a decreasing function of time but its rate of reduction vanishes as 9 approaches 4. All the 
previous treatment of the nonlinearity can be repeated. The only difference in the results is that the timescales td 
and tes are modified. In particular 



td{6) 



1 

2K 



irAk* 9 T(l - 9/4) l 2 D 



4/(4-6) 



(54) 



With this expression one can assess the relevance of the nonlinearity by comparison with the time scale of linear 
theory t l min (9). For 9 < 2 one finds that the nonlinearity is irrelevant [td(&) <C t l min (6)] for 



W 2 k 2 l 2 
VVq k i d 

4a 2 



1 



(2 - 0)r(l - 9/4) 



9 (k* 
1_ 2 \M 



2-9 



r[(2-g)/4] ] 

2T(1 - 6>/4) J 



{4-0)11 



(55) 



The right hand side of the previous inequality is of the order of one for 9 — > 0, in agreement with Eq. (p7|). For 
small 9 it decreases as (h/ko) . For 9 — ► 2 it diverges, as a consequence of the fact that t l min (9) goes to infinity. 
For 2 < < 4 the condition for the irrelevance of the nonlinearity becomes 



W 2 k 2 l 2 
4a 2 



1 



(2 - 0)r(l - 9/4) 



'k* 
2 \hi 



2 9 - 21 



(4-e)/2 



(56) 



The right hand side diverges for 9^2 and vanishes for 9 — > 4, as expected since Tr>(9) diverges in that limit: For 
0^4 the nonlinearity is always relevant. 

A more immediate perception of the meaning of Eqs. ( p5[ ) and ( |56| ) is given by Fig. Jl|, showing the right hand 
sides of the inequalities as a function of 9. Nonlinearity is irrelevant for values of [Wokflr)/ (2a 2 ko)] 2 smaller than 

the function plotted. Except for a small region around 9 = 2, where it diverges (because t l min — > oo), the function is 
always smaller than its value for 9^0 p5[ . This means that correlations increase the effect of the nonlinearity for 
almost all values of 9. 



IV. SOLUTION IN THE NOISY CASE 



So far we have neglected the presence of noise in Eq. (0). We now turn to the study of the problem in presence 
both of deposition and diffusion noise. It will turn out that noise affects the position of the minimum only for small 
initial roughness. For large and intermediate values of Wo the deterministic theory of Sec. [II is sufficient. 

The inclusion of noise in the problem changes Eq. (^) to 



d t S(k, t) = -2 [k^ - f[a(t)]k 2 ] S(k, t) + R{k), 



(57) 



with R{k) = R S + R D k 2 
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The formal solution is 



S(k,t) = S de t(k,t) + S noise (k,t), (58) 
where Sdet(k, t) is again given by Eq. (TL3) while 



Snoise{k,t) = R{k)S det {k,t) ds S- e \(k,s). (59) 



As before the key point is to realize that provided k m /k$ <C 1, after a short transient of duration t$ = l/(2Kfcg) one 
can safely take S de t{k,t) w So exp (— 2nk i t}. Using this expression 

S nolS e(k,t) = [1 - exp {-2nkH)\ . (60) 

With this formula one can compute the additive contributions of the noisy part of the structure factor to a(t) and to 
the roughness W 2 {t) which arise due to shot noise (S) and diffusion noise (D), respectively, 

a(t) = a det (t) + a s (t) + a D (t) (61) 

and 

W 2 (t) = W 2 det (t) + W§ (i) + W% (t) . (62) 

The results are 

a s (t)«^ log (2«**t) (63) 
M*) * (64) 

and 

Wl(t)*nR s] J— (65) 

Wg(t)«^log(2«*St). (66) 

The determination of the temporal evolution of the system is now more complicated than in the noiseless case. There, 
the form of a(t) was always the same and the minimum for W 2 changed depending on the various approximations 
for f[a(t)]. Here, even the expression of a(t) varies in time. Moreover, in the noiseless case, the minimum in W 2 
occurs when the condition (|l^) is fulfilled. Here, there are three independent contributions to the width: A minimum 
may appear much before the instability starts to play any role, simply because of the interplay between the different 
contributions to W 2 . However, since W 2 et is the only decreasing contribution, it is clear that a minimum in W 2 can 
occur only if it already existed in the deterministic case. Noise cannot create a minimum and, as will be shown below 
cannot destroy it, but only shift it to shorter times. We will assume in the following that a minimum in the noiseless 
case exists. 

The complication related to the different contributions to a(t) turns out to be unimportant for physical values 
of the parameters. Inserting the expressions for as and cld in Eq. (0) one sees that, under the physically sensible 
assumptions Id 3> an and fc^ 1 ^> an , /(as) ~ /(«d) ~ a. Hence, even if for long times it may happen that as or ajj 
become larger than a d et, their value is so small that the theory is not affected. 

The analysis of W 2 is instead quite involved. Let us define the timescales 

TRS = 2k^RS = ^ U (67) 
2«W 4 QV 

TRD = ^kjR% = ^T 77 ' (68) 
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TDS= (^J 8^=128^ (69) 

Tfls is the timescale when TVs becomes greater than Wdet', t rd is the timescale when Wd > Wdet', tds is the 
timescale when Ws > Wd- 

Let us carry out the analysis of the interplay of these different timescales assuming, for the moment, that the 
deterministic part is well described by linear theory, i. e. t m i n = t min w 0.18r;. 

For small initial roughness (Q <C 1/2) one finds that trd <C trs <C tds <C t i and this implies that W = Wdet for 
t <C trd, W — Wd for tjid <C t -C r^s, IV = for r^s << <C 77 and IV sa exp(i) for 77 <C £. Hence the minimum 
occurs for t m i n « but before the instability sets in there is another change for t pa i_os- 

If instead Q > 1/2 there are two possibilities. If the initial roughness is not very large (Q 2 y 2 <§C 327r), one has 
tds "C trs "C 77 and IV = W^ e t for i <C trs, W = Ws for trs <C i <C n and IV sa exp(f) for larger times. Otherwise, 
if Q 2 y 2 ^> 327r, then r^s <C 77 <C t^s "C i7j£>: IV = IV^et always, the noise is irrelevant and t min = t l min pa O.I877. 
Notice that the condition Q 2 ?/ 2 32ir is, apart from the numerical factor, the condition (15) of the paper by Krug 
and Rost @. 

We have so far assumed that tf%f n — t l min of the same order of magnitude of 77 . When deterministic nonlinearities 
are strong they increase the value of £^'„, which becomes much greater than 77. As shown above, this starts to happen 
for td 3> T\ which corresponds to <5 2 y 2 3> (8/d/o||) 2 » 1. Then it might in principle happen that f^f n becomes larger 
than trs in the last of the cases above. However, while t^f n grows for large Q as Q (Eq. (0)) or as Q 4 / 3 (Eq. (0)), 
T fls is proportional to Q 2 : r^s always remains larger than t^f n and nothing in the previous discussion changes. 

In Table | is a summary of the different regimes found depending on Q. 



V. NUMERICAL RESULTS 

The results presented above are obtained by analytically estimating the behavior of the large- N equation (Q), which 
in turn is an approximation of the fully nonlinear Eq. (pf). In order to check the validity of these results we have 
solved numerically the fully nonlinear equation for values of the parameters corresponding to the different possible 
regimes of Table [| and compared the results with the numerical integration of the large- N equation pr| and with the 
analytical estimates. The numerical integration was performed by the simple first-order Euler scheme, on a lattice 
of size 512 x 512. The temporal stepsize was chosen to be 1, while the lattice spacing was equal to ir/ko, with fco 
specified in the figure captions. 

We start by considering the limit of very small initial roughness, so that Q = 0.1 and the minimum is due to 
conserved noise roughening the surface (Fig. |^). The minimum width for the fully nonlinear solution occurs for a 
time compatible with the analytical prediction t m i n = trd ~ 1600. Despite having a minimum around the same time 
the solution of the fully nonlinear case and the large-TV approximation differ noticeably for large times. This poor 
agreement is however only apparent and is due to a technical subtlety: The numerical solution of the full equation 
is performed on a square lattice, while the analytical calculations assume a circular Brillouin zone, k < kg. When 
noise is irrelevant, since the structure factor decays exponentially for large wavevectors, the difference in the Brillouin 
zones does not really matter after the initial transient to. In the noise dominated cases instead, the structure factor 
has a power-law tail: The effect of the different Brillouin zones persists in time, cannot be eliminated easily |27| and 
leads to a systematic overestimate of the value of W 2 (t). This is why for long times the numerical solution is not in 
agreement with the large- N result. This problem is most evident in this case dominated by conserved noise, as the 
power-law tail of the S(k,t) is broader. 

In Fig. |^ the value of the roughness of the fully nonlinear case and of the large- N approximation are plotted in the 
case where nonconserved noise dominates. Again, the analytical value t m i n = trs ~ 15700 matches quite well the 
numerical results. 

The same quantities are reported in Fig. || for values of the parameters such that both noise and nonlinearities 
are irrelevant. In this case, the analytical prediction for the position of the minimum is the one provided by linear 
theory t l min pa 0.18r; w 14400. Also in this case, the agreement between numerics and the theoretical prediction is 
good. Notice that for these values of the parameters the naive condition (^6|) for the irrelevance of nonlinearity is 
violated. However, the nonlinear contribution to the current (^) is initially large but rapidly decays, so that it does 
not influence the position of the minimum. The initial deviation from linearity can also be seen in the behavior of 
b(t), plotted in the inset of Fig. ||. 

The position of the minimum is instead determined by the nonlinearity in Fig. |5[ Here linear theory would predict 
t min w 300, while the roughness keeps decreasing until t pa 60000 in reasonable agreement with the analytical 
estimate t m i n « 40000 given by Eq. (j42]). The nonlinear behavior is also evident in the inset, where b(t) is plotted: it 
is proportional to t 3 / 2 as predicted by Eq. (|35|). 
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Finally, in Fig. || the roughness in the most nonlinear case {Q ^> SI d / {a\\y 2 )) is shown. In this case the linear 
theory would predict i TO , n w 0.03, while Eq. ( ^i|) yields t min w 5000. The numerical result t OT j„ w 8000 is again close 
to the estimate provided by nonlinear theory. Consistently b(t) grows as t 2 (Fig. ||, inset). 

In summary, for initial values of the roughness ranging from very small to very large we find that the evolution of 
the fully nonlinear equation is well approximated by the large- N limit and that the analytical estimates found above 
agree with the numerical results. The large- TV limit describes quite precisely this early stage behavior because, up 
to the time when the minimum is reached, the dynamics makes the surface smoother, reducing slope fluctuations 
and making the approximation (Vh) 2 = ((V/i) 2 ) increasingly more accurate. Only after the minimum, when the 
instability takes over, slope fluctuations grow, leading to the breakdown of the large- N approximation. 

VI. DISCUSSION 

In the previous Sections we have carried out a rather complete analysis of the non monotonic behavior of the 
roughness of a surface which evolves according to the continuum equation for unstable growth. It must be stressed 
that an initial decrease of the roughness is not necessarily due to a "rough" substrate, in the sense of surface width 
exceeding some threshold. No matter how small the substrate fluctuations, if they extend to a length scale smaller 
than that of the linear instability, the roughness will initially decrease. Only if substrate fluctuations are limited to 
relatively large length scales, big amplitudes are needed. 

The initial decrease of the roughness is always governed by the relaxational Mullins-like term. The moment when 
this initial decrease ends depends instead crucially on the value of Q (i.e. on the initial roughness): the minimum 
may be accounted for by linear (noisy or deterministic) or nonlinear theory. Interestingly, for realistic values of the 
diffusion length Id and of the ratio y = Ies/Id, the different regimes are nonoverlapping: It is in principle possible 
to see each behavior by simply changing Wq, i. e. the initial roughness. It is clear, of course, that this is not 
necessarily true in practice, since experimentally realizable values of Q are limited. Notice, however, that Q depends 
on two quantities, Wo and ko, and also the variation of the latter could help in expanding the range of variation of 
experimentally realizable values of Q. Moreover the presence of correlations in the initial roughness enhances the 
effect of nonlinearity. 

With regards to the experimental relevance of the present results it is natural to compare them with the recent 
data of Gyure et al. M . The comparison of linear theory with the same data |ll[] pointed out the violation of the 
naive condition for the irrelevance of nonlinearity, Eq. ( p6| ) . As shown above the correct condition for the irrelevance 
is different (Eq. (^)) and turns out to be fulfilled: The substrate is such that nonlinearity does not matter. On the 
other hand, the experimental parameters give Q ~ 110 and y32n/y ~ 450, indicating that nonconserved noise mostly 
dictates the position of the minimum. This conclusion is at odds with the results of Ref . |[l] . The mismatch is due to 
a different treatment of numerical prefactors and should not be taken too seriously: The precision of values determined 
from experimental data is poor and already introduces large uncertainty in the physical parameters. However, even 
if precise experimental data were available a very detailed comparison between theory and experiment would not be 
possible, because formulas linking parameters of the continuum equation (a, k) with physical quantities (Id,Ies) are 
known only in order of magnitude. Improved determination of the numerical prefactors would surely be an important 
contribution to this field of research. 

We have considered here only the form (jj) of the unstable current, which is valid in the limit of small ES barriers 
and does not vanish for finite slopes. Other forms of the current are commonly used for large ES barriers (Ies ^> Id) 
or when the current vanishes for some "magic" value of the slope. The previous analysis can be performed along the 
same lines for these alternative currents. We do not expect qualitatively different results. In particular, no special 
behavior should be induced by the presence of magic slopes, since the initial decay of fluctuations governed by the 
Mullins-like term quickly washes out large slopes independently from the expression of the current. 
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Value of Q 


Position of the minimum 


Relevant effect 


Q< 1/2 
1/2 « Q « V32^/y 
%/32^/y < Q < 8l D /(a ]]y ) 
8l D /(a n y) < Q < 8l D /(a n y 2 ) 
8l D /(a ]} y 2 ) < Q 


TRD = (QV/80^ 

trs = (Q 2 y 2 /327r)r i 
0.18t ; 
0A(tit d ) 1/2 
OMnTDTEs) 1 ' 3 


Conserved noise 
Nonconserved noise 
Linear deterministic 

Nonlinear 

Nonlinear 
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TABLE I. Analytical estimates of the position of the minimum depending on the substrate roughness parameter 
Q — Wo I (a±auko) ■ The other quantities appearing in the Table are determined by the characteristic length scales Id and 
Ies of t he g rowing surface. In particular, y = Ies/Id and r; = An/a 2 w 16F -1 (Id/Ies) 2 ■ For definitions of td and tes see 
Eqs.(|||f. 
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FIG. 1. Plot of the right hand side of equations l55h and (M) vs 6 for h/k = 1/10 and k* /k = 1/100. 




FIG. 2. Double logarithmic plot of W 2 (t) vs t, for a system with k = 1/10, W = 1/100, Q = 1/10, F = 10" 7 , l D = 
and y = 1/100. Here and in all other plots a±_ and an are taken to be equal to 1. For these values of the parameters Q <H 
and conserved noise dictates the position of the minimum. 
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FIG. 3. Double logarithmic plot of W 2 (t) vs t, for a system with k = 1/40, W = 10/4, Q = 100, F = 1, l D = 40 and 
j/ = 1/100. For these values of the parameters 1/2 <C Q <C -\/32ir/y: The minimum is due to nonconserved noise. 
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FIG. 4. Double logarithmic plot of W 2 (t) vs t, for a system with k = 1/100, W = 3, Q = 300, F = 1/50, l D = 265.9 
and jy = 1/10. The inset shows the log-log plot of b(t) along with a line of slope 1. For these values of the parameters 
VS2Ti/y «Q« 8lD/(a,\\y): Linear noiseless theory holds. 
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FIG. 5. Double logarithmic plot of W 2 (t) vs t, for a system with k = 1/100, W = 50000, Q = 5 • 10 6 , F = 10 4 , l D = 10 
and y = 1/1000. The inset shows the log-log plot of b(t) along with a line of slope 3/2. For these values of the parameters 
d>lo/(a\\y) <Q< 8Zu/(o||j/ 2 ) and the nonlinearity is relevant. 
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FIG. 6. Double logarithmic plot of W 2 {t) vs t, for a system with k = 1/1000, W = 10 5 , Q = 10 8 , F = 10 4 , l D = 100 
and y = 1/10. The inset shows the log-log plot of b(t) along with a line of slope 2. For these values of the parameters 
8Z D /(a||y 2 )«Q. 
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